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We determine conditions on the filling of electrons in a crystalline lattice to obtain the equivalent of 
a band insulator - a gapped insulator with neither symmetry breaking nor fractionalized excitations. 
We allow for strong interactions, which precludes a free particle description. Previous approaches 
that extend the Lieb-Schultz-Mattis argument invoked spin conservation in an essential way, and 
cannot be applied to the physically interesting case of spin-orbit coupled systems. Here we introduce 
two approaches, the first an entanglement based scheme, while the second studies the system on 
an appropriate fiat ‘Bieberbach’ manifold to obtain the filling conditions for all 230 space groups. 
These approaches only assume time reversal rather than spin rotation invariance. The results depend 
crucially on whether the crystal symmetry is symmorphic. Our results clarify when one may infer 
the existence of an exotic ground state based on the absence of order, and we point out applications 
to experimentally realized materials. Extensions to new situations involving purely spin models are 
also mentioned. 


I. INTRODUCTION 

Insulating states of matter arise, in clean systems, as 
a result of a commensuration between particle density 
and a crystalline lattice or a magnetic field. Mott in¬ 
sulators are a particularly interesting class, with an odd 
number of electrons in each unit cell. Their low energy 
physics is captured by a spin model with an odd number 
of S' = 1/2 moments in the unit cell. A powerful result 
due to Lieb, Schultz, and Mattis in ID^, later extended 
to higher dimensions by Hastings and Oshikawa^’^, holds 
that if all symmetries remain unbroken, the ground state 
must be ‘exotic’ - such as a Luttinger liquid in ID, or 
a quantum spin liquid in higher dimensions, with frac¬ 
tional ‘spinon’ excitations. These exotic states cannot be 
represented as simple product states, as a consequence 
of long ranged quantum entanglement. This general re¬ 
sult has experimental consequences - indeed no sign of 
magnetic or spatial symmetry breaking is observed down 
to temperatures that are orders of magnitude below the 
intrinsic energy scales in certain materials^, including 
the quasi 2D Mott insulators k-(BEDT-TTF) 2 Cu 2(CN)3, 
/3'Pd(dmit)2 and Herbertsmithite ZnCu3(OH)6Cl2, as 
well as the 3D Mott insulator Na4lr308. Hence if we 
can apply the Hastings-Oshikawa-LSM theorem (col¬ 
lectively referred to as HOLSM) to these systems, a 
strong case is made for an exotic ground state (assuming 
that the effects of disorder can be ignored). However, 
HOLSM invoke spin rotation invariance in an essential 
way, which is typically broken in real materials due to 
spin-orbit coupling. These effects are not small: Her¬ 
bertsmithite has SO(3)-breaking Dzyaloshinskii-Moriya 
terms thought to be on the order of 10% of the Heisen¬ 
berg coupling^’®. In the anti-ferromagnetic hyperkagome 
compound Na4lr3 08, the physics is even dominated by 
spin-orbit coupling effects and charge fluctuation is sig¬ 
nificant^. Physically, the only exact symmetry is time 
reversal (TR) symmetry, and the crystal translations and 


charge conservation which allow us to define the filling. 
Can HOLSM be extended to this physically more realistic 
situation? 

In this work we show that it indeed can, although en¬ 
tirely different theoretical approaches are needed. We 
argue that if a spin-orbit coupled insulator at odd fill¬ 
ing is time-reversal symmetric, its ground state must, in 
a precise sense, be exotic. We introduce two theoreti¬ 
cal approaches that, like the flux threading arguments 
of HOLSM, are non-perturbative, but differs from them 
in that conservation of spin is not assumed. The first is 
an entanglement based approach that allows us to prove 
that symmetric, gapped and unfractionalized insulators - 
the interacting analog of a band insulator, with a unique 
ground state on both the plane and torus - are only al¬ 
lowed at even integer fillings v = 2m. For brevity we re¬ 
fer to such symmetric short-range entangled insulators as 
‘sym-SRF’ insulators. A corollary of this result is that at 
odd integer fillings, Mott insulating phases must either 
break a symmetry or have a ground state degeneracy 
on certain geometries due to other, more exotic, mech¬ 
anisms. A special case of this result in ID spin models 
was previously discussed in Ref. 7. Here we will extend 
it to higher dimensions and allow for charge fluctuations. 

This constraint on filling arises even when translations 
are the only spatial symmetries. What if additional sym¬ 
metries are present, such as the 230 space groups of 3D 
crystals? It turns out that additional constraints ap¬ 
pear only for the non-symmorphic space groups, where 
the minimal filling at which a sym-SRF insulator arises 
is at least v = 4. We find lower bounds on the mini¬ 
mal filling for all 157 non-symmorphic space groups, and 
these bounds are shown to be the tightest possible for 
a large majority of them. Earlier results on noninter¬ 
acting band structures®“^° had pointed out that in non- 
symmorphic crystals there are required band touchings 
leading to larger minimal fillings. In Refs. 11 and 12 
this was generalized to interacting systems using flux 
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threading arguments. However, similar to the HOLSM 
arguments, these need to assume spin rotation invari¬ 
ance, and typically do not lead to useful constraints in 
their absence. Our entanglement based argument allows 
both for strong interactions and broken spin rotation in¬ 
variance, while preserving time reversal symmetry. The 
best bounds are obtained using a different approach, by 
defining the system on a closed manifold that is locally 
flat and hence is locally indistinguishable from Euclidean 
space (the Bieberbach spaces). These are generalizations 
of the torus - where instead of identifying points related 
by translations, one uses non-symmorphic elements, say 
glides or screws, to obtain a compact manifold. By expos¬ 
ing an obstruction for systems defined on such spaces, we 
show that the sym-SRE phase, which should be sensitive 
only to local physics, cannot appear at certain fillings. 

Finally we show there are special cases where obstruc¬ 
tions to a sym-SRE insulator occurs despite being on 
symmorphic lattices at even integer filling. These only 
occur in spin models without charge fluctuations, where 
reflection symmetry relates two entanglement cuts that 
can enclose an odd number of sites. 


OVERVIEW 

Before presenting our detailed arguments, we first give 
an intuitive description of the key ideas and results. 

In this work we prove a general constraint on the min¬ 
imum number of electrons per primitive unit cell 
required to realize a sym-SRE phase: 

Consider a system of electrons with Kramers degener¬ 
acy. We assume the U(l) (particle number) symmetry, 
a space group symmetry (including lattice translations), 
and the time-reversal symmetry. Then, to obtain a sym- 
SRE phase the average number of electrons per primitive 
unit cell (i.e. filling) has to be an even integer no less 
than two: Vmin > 2. Furthermore, if the space group is 
non-symmorphic, Vmin must be an even integer no less 
than four. 

To warm up, we first ignore charge fluctuations and 
consider an infinite, periodic ID chain with one localized 
spin-1/2 moment per cell. For an conserving chain at 
half filling (including SO(3) invariant chains), the LSM 
theorem asserts that the ground state must be either gap¬ 
less or break a symmetry. A more general version of LSM 
is known in ID^: when each unit cell carries a projective 
representation of a symmetry of the system (or Kramers 
degeneracy), a sym-SRE phase is forbidden. 

To illustrate the ideas discussed in Ref. 7, consider a 
spin-1/2 chain with TR symmetry T. Suppose the sys¬ 
tem is in a sym-SRE phase, and consider an entangle¬ 
ment cut at some bond x, which determines a Schmidt 
decomposition of the ground state. Due to the presence 
of a gap, the high-weight Schmidt states have a discrete 
spectrum,which we label by a. The Schmidt states 
\ot)L are all either Kramers singlets, with T^ja)/, = |a)i, 
or doublets, with = —1. This is because in the ab- 
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FIG. 1. Schematic Schmidt decomposition of an infinite 
periodic chain with one spin-1/2 per cell and time-reversal 
symmetry, (a) Two entanglement cuts, related by a lattice 
translation, enclose one spin-1/2 moment. If the system is 
in a sym-SRE phase, the Schmidt states have the same 
eigenvalue at the same cut and has to change from, say, = 
-1-1 (blue) at cut x to = — 1 (red) at cut x -I- 1, which 
is inconsistent with translation symmetry. The system can 
circumvent such obstruction by either becoming gapless or 
spontaneously breaking symmetry (such that the symmetric 
state is a cat state), (b) Pictorial description of a VBS cat 
state. 


sence of long-range order, all Schmidt states differ by the 
action of operators localized near the cut, and local op¬ 
erators never switch the value of T^. Now advance the 
entanglement cut by a lattice translation to bond a: -I- 1, 
as illustrated in Fig. 1(a). On the one hand, we have 
added one more spin-1/2 to the Schmidt states, so their 
eigenvalue must change by —1; on the other hand, 
the two cuts should be equivalent due to translation in¬ 
variance, so their eigenvalue should not change. We 
thus arrive at a contradiction: the Kramers degeneracy 
carried by each unit cell translates into an obstruction to 
a sym-SRE phase. 

In ID there are two ways to circumvent the obstruc¬ 
tion. If the system is gapless, the Schmidt weights form a 
degenerate continuum and it is no longer well-defined to 
track the eigenvalue of any particular Schmidt state. 
If the system spontaneously breaks a discrete symmetry, 
one can still form a symmetric combination of the ground 
states, but it will be a Schrodinger cat state with long 
range order. For a cat state the Schmidt states look dif¬ 
ferent from each other arbitrarily far away from the cut, 
and the assignment of eigenvalues may differ across 
different Schmidt states at the same cut. For example, 
for the valence-bond solid (VBS) state in Fig. 1(b), we 
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assign = — 1 if the entanglement cut dissects a sin¬ 
glet pair, and = +1 otherwise. Consequently the 
eigenvalue depends on the Schmidt state and our argu¬ 
ment breaks down. 

To obtain similar bounds in higher-dimensional peri¬ 
odic systems, one can extend the above argument by 
suitably imposing periodic boundary conditions with odd 
circumference in all but one dimension to form a quasi- 
ID cylindrical geometry, then apply the ID argument. 
In dimensions higher than ID, if the system is gapped 
the resulting ground state degeneracy dX v < can 
either be due to spontaneous symmetry breaking, or 
due to topological order. For topological orders with 

< ^'mini the topologically degenerate ground states as¬ 
sociated with periodic boundary conditions may be re¬ 
lated to each other by the symmetries, which is a subtle 
form of symmetry-breaking. For instance, this implies a 
TR-symmetric system of a Kagome lattice, having three 
spin-1/2 moments per cell, cannot be in the sym-SRE 
phase. 

HOLSM is only one example of a much larger class of 
constraints we obtain using this program: a sym-SRE 
phase is actually prohibited whenever two entanglement 
cuts are related by some spatial symmetry (such as re¬ 
flection, glide, or screw) and they enclose a projective 
representation. In particular, tighter bounds can be ob¬ 
tained when the enclosed volume is only a fraction of the 
unit cell. Later in this work we give a concrete model of 
localized spin moments in which the two cuts are related 
by a reflection. Consequentially a sym-SRE phase can 
be forbidden even when there are two spin-I/2 moments 
per unit cell. 

In real materials, however, the spin-1/2 moments typ¬ 
ically come from spinful electrons that are not strictly 
localized and spin-orbit coupling can be significant. 
Naively, one might think any lower bound on the required 
number of localized spin-1/2 moments simply morphs 
into one on electron filling. Surprisingly, this is not 
true: there are localized spin models, like the reflection- 
symmetric one mentioned above, for which a sym-SRE 
phase ceases to be forbidden once charge fluctuation is 
allowed. 

To understand real materials, it is therefore impor¬ 
tant to extend HOLSM-type theorems to electronic sys¬ 
tems with both charge fluctuation and spin-orbit cou¬ 
pling. We address this problem in this work assuming 
only particle number conservation and time-reversal sym¬ 
metry. Our results indicate that an electronic system 
can be in a sym-SRE phase only when the filling is an 
even integer no less than two. For most non-symmorphic 
space groups, the presence of fractional translations in 
the glides or screws allows us to obtain an even tighter 
bound of t'niin > 4. Two non-symmorphic space groups, 
however, are exceptional and do not contain any intrin¬ 
sically non-symmorphic elements^^; they instead feature 
a ‘Borromean ring’ of screws, and our entanglement ar¬ 
gument does not distinguish them from the symmorphic 
ones. 


TABLE I. Summary of minimal filling i/min for sym-SRE 
insulators for ‘elementary’ space groups. 




Minimal Filling 

(Conway) 

ITC No. 

Sym. " 

A.I.*’ 

Ent.'’ 

Bbb.'^ 

Manifold name^® 

1 ^ 

(translation) 

2 

2 

2 

Torus T® 

4 

2i 

4 

4 

4 

Dicosm 

144/145 

3i/32 

6 

6 

6 

Tricosm 

76/78 

4i/43 

8 

8 

8 

Tetracosm 

77 

42 

4 

4 

4 


80 

4i-bBCTf 

4 

4 

4 


i69/i70 

61/65 

12 

12 

12 

Hexacosm 

171/172 

62/64 

6 

6 

6 


173 

63 

4 

4 

4 


19 

2i2i 

8 

4 

8 

Didicosm 

24 

2i2i-bBCT 

4 

2 

4 


7 

glide 

4 

4 

4 

U* Amphicosm 

9 

glide 

4 

4 

4 

2"'* Amphicosm 

29 

glide-l-2i 

8 

4 

8 

U* Amphidicosm 

33 

glide-l-2i 

8 

4 

8 

2"'* Amphidicosm 


‘Sym.’ lists the key symmetry elements of the space groups 
that are relevant to our discussion. 

The minimal electron filling required to form a TR-symmetric 
atomic insulator 

“ Results in section ‘Extension to 3D symmorphic and 
non-symmorphic crystals’ 

Results in section ‘Alternative method: Constraining the filling 
by defining sym-SRE insulators on compact flat manifolds’; for 
4 of the elementary NS space groups (No. 19, 24, 29 & 33) this 
gives a tighter bound 

® For completeness we include results for the ‘elementary’ 
symmorphic space group 

^ BCT; Body-centered translation; the minimal filling is therefore 
half of the corresponding primitive space group 

To handle them, we appeal to the following physical 
intuition: a sym-SRE phase is sensitive only to local 
physics, so should have a unique ground state on any 
large lattice geometry that looks locally identical to the 
translationally invariant infinite plane. These geometries 
are the lattice analogs of compact flat 3-manifolds, and by 
systematically studying all the compact flat 3-manifolds 
we obtain the tightest possible lower bounds for a set of 
‘elementary’ space groups that we list in Table 1. Using 
the group-subgroup relations, we manage to derive lower 
bounds on the sym-SRE electron filling for each of all 230 
space groups (Fig. 2 of Appendix). These lower bounds 
are provably tight for all but ten space groups. 


A LSM-type theorem for interacting spin-orbit 
coupled chains 

Here we prove that a TR-invariant ID interacting SOC 
chain at odd filling ^ is either gapless or breaks symmetry 
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(translation or TR). The argument presented here para¬ 
phrases a more precise proof using the language of matrix 
product states in Appendix. We proceed by assuming the 
ground state |4') is symmetric and short-range correlated, 
and derive a constraint on the filling v. Let x label unit 
cells and Ti be the translation x —>■ x -\- 1 . Each unit 
cell may contain many sites and orbitals. We Schmidt 
decompose |4') across a bond x S (xq — l,xo), 

1 ^) = ( 1 ) 

a 

with Schmidt weights s^q. Assuming the phase is 
gapped, obeys an area law, so the spectrum s^a (and 
hence the label a) is discrete. Each Schmidt state is an 
eigenstate of the number operator, but it is subtle to dis¬ 
cuss their eigenvalues (the total charge to the left) as they 
are ill-defined in the thermodynamic limit. Instead, the 
well-defined quantity is the charge relative to the mean 
filling, Qx = Nx — v. We have 

^ ' Qx\ci) X L ~ Qo‘\^) X L ^ Qa = 'y ^ {^\Qx\oi)x L ■ 

x<,x x<.x 

( 2 ) 

The sum is well conditioned because {a\Qx\a) ^ ^ ^ Q as 
X ^ —oo exponentially quickly for a gapped state. This 
follows from the fact that the Schmidt states differ from 
the ground state only in the vicinity of the cut and the 
ground state is assumed not to have charge-density wave 
(CDW) order. 

Translation invariance Ti |'I') = implies that the 
Schmidt decompositions across different bonds are re¬ 
lated by Sxa = Sx+la, Tl W) x L / R = \^) x+1 L / R^ 
two translation-related states having the same charge ex¬ 
cess Qa- The Schmidt states at the two cuts are related 
by the addition of the intervening sites, spanned by states 

\p)x- 

\<^)x+lL='^Kfl\p)x\l^)xL- (3) 

Choosing a charge eigenbasis for the sites, Qx \p)x = 
Qp \p)^ with Qp S Z — z/, charge conservation implies that 
whenever ^ 0 we have Qa = Qp + Qp- However, so 
long as the state is short-range correlated, Qa — Qp GZ, 
since any two Schmidt states differ by the addition or re- 
arrangment of some particles near the entanglement cut. 
This is a contradiction unless ly G Z (the HOLSM theo¬ 
rem). If 14/) is a cat state of two CDWs, so is symmetric 
but not short-range correlated, the differences Qa — Qp 
no longer need to be an integer and our argument does 
not apply. 

When TR symmetry is incorporated, each Schmidt 
state is either part of a Kramers singlet or doublet: 

= (r^)« = ±i- (4) 

Of course {T^)a is related to (—1)‘5“, but we have to 
recall that is dehned to be (—1)'^“’, not (—1)'^“. The 


divergence of the total charge results in an ambiguous 
but (assuming short-range correlations) a-independent 
phase: 

(r2)a=e*^(-l)0“. (5) 

Again, because of translation invariance, the charge ex¬ 
cess Qa and the TR character {T‘^)a are independent of 

X- 

Now consider the action of on the decomposition 
(3): 

= Y.Blp{f^ b).o)(r" \P)xl) 

= \p).omxL 

p,p 

= ( 6 ) 

where we used Qa = Qp + Qp and Eq. (3) again in the 
last step. Note the important factor of (—I)'^ compared 
to (T^)a in Eq. (5). It follows that (—1)"^ = 1 and hence 
1 / G 2Z. 

\iv ^ 2Z, one of our assumptions must be violated: ei¬ 
ther the phase is gapless, the state breaks a symmetry, or 
it is symmetric but a long-range correlated cat state. For 
example, for a gapless system the entanglement Hamil¬ 
tonian is also gapless: all Sa -G 0, the Schmidt states 
form a continuum, and the symmetry properties of the 
Schmidt states are ill-defined. Alternatively, for a VBS 
cat state in Fig. 1(b), the ambiguous phase is no 
longer a-independent, but instead will depend on which 
of the two VBS sectors the Schmidt state belongs to. 

A. Extension to 3D symmorphic and 

non-symmorphic crystals 

The above argument can be extended to 3D systems^® 
using the quasi-lD (generalized) cylinder K x T^. Given 
a 3D crystal with the primitive lattice vector ai, 2 , 3 ) we 
introduce a periodic boundary condition that identifies x 
with X -G N 2 a 2 and x + N^as- We are then left with a 
single ‘long’ direction with infinite extent along ai and 
a translation symmetry Ti. Formally we can view the 
cylinder as a ID-chain with a ID-unit cell that comprises 
N 2 N 3 of the 3D unit cells. Therefore, choosing N 2 , A 3 to 
be odd, the constraint on ID hlling z^id = {N 2 N 3 )iy G 2Z 
implies v g2Z- 

The extension implicitly assumed the following: if the 
phase is trivial, for sufficiently large N 2 ,N^ the resulting 
cylinder must also have a unique ground state. Intu¬ 
itively, a sym-SRE phase cannot detect the global topol¬ 
ogy - we will comment further on this later. Similar as¬ 
sumptions were invoked in Oshikawa’s argument^. 

The lower bound Vmin = 2 is in fact the tightest bound 
for all 73 symmorphic space groups in 3D (unless extra 
constraints, like specifying a lattice realization, are im¬ 
posed). A space group is said to be symmorphic if every 
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symmetry element can be decomposed into a product of 
a lattice translation riidi (jii S Z) and a point group el¬ 
ement. Such a decomposition implies that putting a site 
at the ‘origin’ of the point group within each unit cell 
generates a lattice that respects all the symmetries. For 
symmorphic space groups, therefore, one can trivially re¬ 
alize an atomic insulator at = 2 by putting a pair of 
electrons into a singlet state at the highest symmetry 
point in each unit cell. Since we can explicitly construct 
a sym-SRE state aX i/ = 2, the bound iZmin = 2 is tight. 
Note that such a high-symmetry point is a special in¬ 
stance of the allowed positions of sites compatible with 
the space group, known as the ‘Wyckoff positions’. A 
set of sites in a Wyckoff position (labeled by a ‘Wyckoff 
letter’) transforms within itself under all symmetry ele¬ 
ments, and can be thought of as the admissible positions 
of the atoms. 

The situation is very different when the space group 
is non-symmorphic (NS), for which given any fixed ori¬ 
gin at least one symmetry element contains a fractional 
translation. As a direct consequence, each unit cell must 
contain more than one site. Since an atomic insulator is 
formed by having a singlet pair of electrons per site, the 
upper bound on v^nin per cell can be inferred from the 
Wyckoff position with fewest sites. This is listed in the 
‘A.I.’ column of Table I. 

To systematically study all 157 NS space groups, we 
focus on those 18 space groups listed in Table I, which 
we refer to as ‘elementary’ in the following. All other 
NS space groups contain at least one of these elemen¬ 
tary groups as a subgroup. (More precisely, here we are 
referring to so-called translationengleiche subgroups or t- 
subgroups^^, which retain the primitive unit cell of the 
entire group.) Thus, a lower bound on the elementary 
space groups immediately provides a lower bound for 
the others, though further work is required to determine 
whether the inferred bounds are tight. 

All elementary NS space groups, except for No. 24, are 
made NS by a n^-screw (27r/n-rotation Cn followed by 
m/n-translation) or a glide (a mirror reflection followed 
by a half translation), which we denote by G and list in 
the ‘Sym. elements’ columm of Table I. No. 24 is known 
as an ‘exceptional’ NS space group and will be addressed 
later. 

To establish a dimensional reduction as before, we have 
to specify our choice of the primitive lattice vectors 01,2,3- 
We require the following: (i) the NS operation G is rep¬ 
resented as G = g-d‘nT'/n)P-aiJ^^ Here, X is the rotation 
Gn about an axis parallel to di for Um-screws, or the 
mirror about a plane that contains Oi for a glide (where 
m/n = 1/2). (ii) The plane spanned by 02 and 03 is 
invariant under X. For these elementary NS groups, one 
can check case by case that such a choice is possible. 
(For No. 80, one has to use the conventional cell, which 
contains two primitive unit cells.) 

As before, we take a cylinder geometry K x by in¬ 
troducing the periodic boundary condition for A^ 2 a 2 and 



FIG. 2. Entanglement cuts enclosing fractions of the unit cell, 
(a) In the quasi-lD setup, two entanglement cuts related by G, 
say a 2i screw, enclose a half-integral volume, (b) Distorted 
checkerboard lattice. Two entanglement cuts are related by 
the reflection symmetry. 

N^ds- We then replace the translation Ti in the above 
argument with the NS operation G. The action of X 
is purely ‘onsite’ in the ID picture and has no effect on 
our ID argument, since X merely rotates or reflects the 
plane spanned by 02 and 03 , as illustrated in Fig. 2(a). 
The volume enclosed by the two cuts related by G is m/n 
times smaller than the symmorphic case, so we now have 
Pm = (m/n)A 2 A 3 Z/. Requiring G 2Z, we prove the 
Pmin listed under column ‘Ent.’ in Table I. (Remember 
that the conventional cell of No. 80 contains two prim¬ 
itive unit cells, so we have to divide the naive bound 8 
by two to get the filling bound of 4.) For most of the 
elementary NS space groups, the lower bound obtained 
here coincides with the upper bound from the atomic 
insulator limit, and hence is also the tightest bound. 

Curiously, the argument here does not produce a 
tighter bound for No. 24. The difficulty is that No. 24 
contains three orthogonal 2 i screws, and if any one is ab¬ 
sent the crystal is symmorphic. Therefore, we need a new 
machinery that makes use of multiple NS operations at 
once. Similarly, the bounds obtained for No. 19, 29 and 
33 are not the tightest because they are also generated 
by multiple NS operations. 


Alternative method: Constraining the filling by 
defining sym-SRE insulators on compact flat 
manifolds 

The above argument based on the ID Schmidt decom¬ 
position can only utilize a single glide or screw, whereas 
multiple such non-symmorphic operations can be present 
in a space group. Here we present a different line of ar¬ 
gument to overcome this challenge. 

We employ the following physical principle: if a Hamil- 
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tonian defined on an infinite lattice A with space group 
SG is in a sym-SRE phase, it will have a unique ground 
state on any compact lattice geometry in which a) the 
Hamiltonian is locally indistinguishable from that on A, 
b) the particle density is that of the thermodynamic 
limit, and c) the width of all compact directions is suffi¬ 
ciently larger than some scale comparable to the correla¬ 
tion length. The lattice geometries we have in mind are 
generalizations of the torus; in the continuum we would 
say they are ‘fiat’. Intuitively, a sym-SRE state is only 
sensitive to local physics, so cannot detect the global ge¬ 
ometry so long as each patch is locally identical to A. 

Let us first use this principle to reproduce our con¬ 
straint on V for symmorphic crystals. To that end, it 
is sufficient to put the system on a torus with odd 
circumferences , TVs in every direction, which is a 

well defined procedure since the Hamiltonian is transla- 
tionally invariant. If the filling v is odd, the total number 
of electrons A’e = is also odd and Kramers de¬ 

generacy is unavoidable. Yet according to the principle, 
a sym-SRE phase should have a unique ground state: it 
is short-range correlated, so cannot detect the parity of 
Ni once Hence, to realize a sym-SRE phase, v 

must be even. 

The advantage of this argument is that one may get 
a stronger constraint on v by using manifolds other that 
the torus. To this end, we discuss how to define a lattice 
model on a fiat manifold obtained by ‘modding-out’ 
by a subgroup T of the SG. Let us first formalize the 
familiar example of putting a translation invariant ID 
system onto a ring S^. If the circumference of the ring is 
N, we identify sites x with x + nN (n € Z) and operators 
Ox with Ox+nN- Denoting T as the group generated 
by an A^-step translation, one can formally rewrite this 
identification rule (~r) as 

X ~r gix), Of ~r gOfg-^ = (7) 

for g S r. Here, g is the quantum operator corresponding 
to 5 , O J runs over all the operators at x, and Ug is a 
unitary matrix representation which could account for 
any action of g on the internal degrees of freedom, as will 
eventually be required for spin-orbit coupling. Eor the 
usual notion of translation, Ug = 1. The Hamiltonian 
for the ring is obtained from that of the infinite chain by 
identifying terms using the equivalence ~r- 

The same procedure applies to more complicated space 
group symmetries: given a 3D lattice Hamiltonian sym¬ 
metric under T C SG, we impose the equivalence relation 
~r in order to obtain a lattice model on the space K^/T. 
There are two requirements on T in order to obtain a well- 
defined lattice Hamiltonian. Eirst, T must be fixed-point 
free (also called ’Bieberbach’), as otherwise the manifold 
will contain conical singularities or edges, so won’t be lo¬ 
cally indistinguishable from the plane (see Appendix for 
more details). Second, since Ug also encodes the non¬ 
trivial action of glides and screws on the internal spin, T 
has to be compatible with certain consistency conditions 



FIG. 3. Illustration of a Didicosm obtained by modding out 
by a suitably chosen subgroup of space group No. 24. This 
space group is generated by two orthogonal screws. The black 
box represents a fundamental domain, and the blue and red 
boxes are the image of the fundamental domain under the two 
screws. The letters P, Q and R and their orientation indicate 
pairs of faces that are identified. 


(namely the existence of the lattice analog of Spin and 
Pin“ structure), which we will return to shortly. 

There are 10 compact flat manifolds (also known as 
Bieberbach manifolds) in 3D, and each of them can be 
obtained by modding-out by a Bieberbach group T^^. 
Eor example, when T is generated by a glide and two 
other orthogonal translations we obtain a manifold dif- 
feomorphic to (the Klein bottle) x , the so-called ‘first 
amphicosm.’^® When T is non-symmorphic, the number 
of unit cells contained in the resulting manifold K^/T can 
be fractional, and will lead to a tighter bound. 

For example, to derive the constraint for space group 
No. 24, we put the system on a manifold obtained from 
a subgroup T spanned by two suitably chosen screws. 
The resulting manifold (Fig. 3) can contain a half-integer 
number of primitive unit cells of No. 24. Specifically, the 
number of electrons will be Ne = vNiN 2 {‘ 2 .N^ + l)/2 for 
Ni G Z. Choosing N 1 N 2 to be odd, we must have 1 / G AZ 
to avoid a Kramers degeneracy. 

By engineering T for each of the elementary space 
groups, we have obtained the tight bounds r'min, as sum¬ 
marized in Table 1. Furthermore, using group-subgroup 
relations described in detail in Appendix, we determine 
the set of all possible fillings that are allowed to realize a 
sym-SRE phase for each of all 157 NS space groups, which 
we include in Appendix. Except for 10 NS space groups 
(No. 73, 106, no, 133, 135, 142, 206, 220, 228 and 230), 
these fillings are necessary and sufficient - there is indeed 
a noninteracting band insulator at these filling^®. For the 
10 NS groups mentioned, however, our method here indi¬ 
cates a sym-SRE insulator is possible at ^ = 4, whereas 
one can in fact show a band insulator at the same filling 
is impossible^®. Therefore, there are two possibilities — 
either one can actually show a stronger constraint that 
prohibits this filling (by utilizing point group symmetries 
of the space group, which we did not use in this work) or 
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there indeed exists an interacting example of sym-SRE 
phase at this filling. We leave this open problem to future 
work. 

In closing we note the following technical point. It is 
well known that consistently defining spinors in the con¬ 
tinuum requires the manifold to admit ‘Spin’ or ‘Pin^’ 
structure, so we should investigate whether this con¬ 
straint arises in the SOC lattice setting (Pin^ structures 
extend Spin structures to non-orientable manifolds; the 
± sign indicates that the square of a reflection is ±1). 
The subtlety arises because SOC spinful fermion opera¬ 
tors fscr transform under a double group of the SG [e.g., 
SU(2) rather than SO(3) for the proper rotation part 
of the SG], leading to sign ambiguities in the choice of 
Ug in Eq. (7). If the signs are not fixed properly, the 
Ug will generate a projective representation of P, i.e., 
^gi ^92 = ■ When the Ug are projective, the equiv¬ 

alence relation is ill-defined; we do not know if we 
should identify Ug^g^fg^g^^g) with or -fg. There¬ 
fore, we have to carefully choose the sign of each Ug in 
such a way that Ug forms a linear (nonprojective) repre¬ 
sentation of r. This choice may not be unique, and as 
detailed in Appendix each choice corresponds to a dif¬ 
ferent Spin / Pin“ structure^®. Fortunately, all 3D flat 
manifolds admit either a Spin or Pin“ structure, so the 
proper choice of Ug is always possible. 

Localized spin models: New symmetry constraints 

So far we have focused on universal constraints on ly 
with the input of space groups only. We now show that, 
given specific positions of sites in the unit cell, there is 
an entirely new class of constraints. Unlike the previous 
section where we could utilize only fixed-point free sym¬ 
metry operations, here we use a point-group symmetry 
to probe the location of the sites in the unit cell. For this 
result we must restrict to the strongly Mott-insulating 
regime where the system reduces to S' = 1/2 moments at 
lattices sites. In the presence of certain reflection sym¬ 
metries, even for symmorphic crystal lattices the mini¬ 
mal filling for a trivial paramagnet can be r'min = 4, not 
i^min = 2. In particular: If an S = 1/2 magnet has an 
odd number of moments per unit cell lying on a reflection 
plane, the system is either gapless, or the ground state is 
not-unique. The ‘tighest’ bound is shattered because here 
we assume a more constraining scenario, namely we as¬ 
sume specific lattice realizations of the space group sym¬ 
metries where the electrons are strongly localized. 

As a concrete example we consider the distorted 
checkerboard lattice shown in Fig. 2(b) with S' = 1/2 mo¬ 
ments localized on every lattice site. This is an example 
of a lattice with the symmorphic wallpaper group symme¬ 
try No. 3 {pmi), which contains a mirror {x, y) — {—x, y). 
There are two inequivalent mirror planes a: = 0 and 
X = 1/2, and we consider a lattice formed by putting 
one site on each of the two planes in each unit cell as 
shown in Fig. 2(b). No symmetry relates the moments 


on the two sub-lattices a and b (corresponding to Wyckoff 
letters a and b in Ref. 17), and the crystal is symmorphic, 
so our previous results do not imply a constraint. 

Now we consider putting the system onto a cylinder 
with an odd circumference in the y direction. Consider 
an entanglement cut at anywhere between x = 0 and 
1/2. Under mirror reflection about the plane x = 1/2, 
the two entanglement cuts enclose an odd number of sites 
in sublattice b and hence an odd number of S' = 1/2 mo¬ 
ments. Now suppose that each S = 1/2 moment carries 
a projective representation of some symmetry, say TR or 
SO(3). Then these symmetry-related entanglement cuts 
yield corresponding Schmidt states carrying different rep¬ 
resentations. This rules out the possibility of having a 
trivial ground state at filling v = 2. 

Note, however, if charge fluctuations are allowed, then 
we can localize both electrons to one or the other lattice 
site and thereby obtain a sym-SRE insulator. Therefore 
it was important that we consider the pure spin model. 
Another route to preserving the obstruction is to consider 
Hamiltonians with particle-hole symmetry, obtained for 
instance, by only allowing hopping between bipartite lat¬ 
tice sites. This also effectively constrains the filling on 
sites. 


APPLICATIONS AND CONCLUSION 

Let us discuss some physical applications of our result. 
First, consider the statement - a Mott insulator with an 
odd number of electrons per unit cell must either be gap¬ 
less or, if gapped, must display a ground state degener¬ 
acy. This is true even when only time reversal symmetry 
is present and spin rotation invariance is absent. The 
Mott insulator has a net Kramers doublet per unit cell, 
which is the analog of having a half odd integer spin per 
unit cell. This is of practical relevance in real materi¬ 
als, where spin-orbit coupling always breaks the full spin 
rotation symmetry while time reversal continues to be a 
physical symmetry. Let us consider the possible ground 
states of a Mott insulator with just time reversal sym¬ 
metry in D = 2. The conventional ground states either 
spontaneously break time reversal symmetry, or spatial 
translation symmetry. In both cases ground state de¬ 
generacies result due to symmetry breaking. However, 
any other cause for ground state degeneracy, or gapless¬ 
ness, must have an ‘exotic’ origin. For example, ground 
state degeneracy that is not attributable to symmetry 
breaking implies there is no local operator that could 
distinguish the different states. This is the hallmark of 
topological order, as in a quantum spin liquid, where 
emergent excitations with self or mutual anyonic statis¬ 
tics are expected. Similarly, the conventional origin of 
gapless excitations in a phase, either a Fermi sea of elec¬ 
trons or Goldstone modes are both ruled out since we are 
dealing with a symmetric insulator. The materials Her- 
bertsmithite (ZnCu 3 (OH) 6 Cl 2 ) is a Mott insulator with 
spin 1/2 moments on a Kagome lattice^^. The ground 


state is believed to be symmetric since no thermody¬ 
namic phase transition occurs on lowering the temper¬ 
ature to a small fraction of the exchange energy. In ad¬ 
dition to Heisenberg exchange terms, spin-orbit coupling 
induced Dyaloshinskii-Moriya interactions lower the spin 
rotation symmetry®. Hence one cannot directly apply the 
Oshikawa-Hastings type argument to this system. How¬ 
ever, since time reversal is preserved, we can use the re¬ 
sults of the present paper to conclude that the resulting 
ground state must be an unconventional phase of matter 
such as a gapped quantum spin liquid. 

A more surprising application of our results concerns 
the truly three-dimensional material Na 4 lr 308 , in which 
the iridium atoms form a hyperkagome lattice with 12 
sites per primitive unit celP^. Similar to Herbertshimite, 
no thermodynamic phase transition is observed when 
approaching the zero-temperature limit and the ground 
state is believed to be symmetric. In a simplified picture, 
the system can be described as an anti-ferromagnetic 
Heisenberg model with one spin-1/2 moment on each of 
the iridium sites. In reality, however, both spin-orbit cou¬ 
pling and charge fluctuation are significant in Na 4 lr 308 ^. 
This precludes any application of the original HOLSM 
theorem. While our results have taken into account 
such generalizations, the electron filling here is even and 
so the argument above, utilizing the odd electron fill¬ 
ing in Herberstmithite, is also inapplicable. Nonethe¬ 
less, we can still rule out a sym-SRE insulating phase by 
incorporating the non-symmorphic space group symme¬ 
tries: Na4lr308 crystallizes in the space group No. 213 
(P4i32)^^, which admits a sym-SRE insulator only when 
iz S 8Z (Appendix Table I) but v = 12. 

One caveat though is the presence of impurities strictly 
speaking, destroys translation invariance which is crucial 
to defining filling^^. In addition, for some systems certain 
atoms are not perfectly ordered and the space group sym¬ 
metries are, strictly speaking, only approximate. This 
is the case in the hyperkagome Na 4 lr 308 , in which the 
sodium atoms sit at Wyckoff positions 4a and 12d with 
only 75% occupation^^. Hence all statements are about 
a putative clean system in which the relevant degrees of 
freedom are symmetric under the space group considered. 

As a second example consider the Dirac semimetal de¬ 
scribed in Refs. 23 and 24. A simple model of this state 
is the Fu-Kane-Mele^® model of electrons on a diamond 


lattice at half filling, which leads to Dirac cones at face 
centers of the Brillouin zone. In a noninteracting electron 
picture, the stability of this interesting electronic struc¬ 
ture arises from the non-symmorphic space group of the 
diamond lattice, which leads to a four dimensional irre¬ 
ducible representation which forms the apex of the Dirac 
cone. Can strong interactions gap these nodes and lead 
to a sym-SRE insulator? Since the diamond lattice has 
two sites in the unit cell, the filling of two electrons is 
compatible with a sym-SRE state in a symmorphic crys¬ 
tal. However, given the non-symmorphic diamond lat¬ 
tice, the minimum trivial filling is four - implying that if 
the Dirac nodes are gapped the system necessarily breaks 
symmetry or is topologically ordered. 

In conclusion, we have obtained the allowed sequences 
of filling for which SOC-coupled interacting electrons 
can form a sym-SRE insulator for all elementary space 
groups. For symmorphic lattices this happens to be 
an even number of electrons per unit cell, while for 
non-symmorphic lattices the sequence is more restric¬ 
tive and demands at least four electrons per unit cell. 
The minimal fillings for a variety of non-symmorphic 
space groups were derived, however a complete explo¬ 
ration of the tightness of these lower bounds is left for 
future work. Similarly, the exploration of constraints on 
topological order and symmetry fractionalization for non- 
symmorphic lattices, analogous to the analysis for Mott 
insulators^®, is an interesting open problem. 


ACKNOWLEDGMENTS 

MZ is indebted to conversations with Mike Freedman, 
Matt Hastings, Mike Hermele, and Meng Cheng and 
thanks Sid Parameswaran for the introduction to NS 
space groups. AV thanks S. Parameswan, Ari Turner and 
Dan Arovas for an earlier collaboration on related top¬ 
ics, Leon Balents for insightful discussions and support 
from NSF DMR 1206728 and the Templeton Founda¬ 
tion. HW thanks Takahiro Morimoto, Yohei Fuji, Masaki 
Oshikawa, and Hitoshi Murayama for useful comments. 
HCP was supported by a Heilman graduate fellowship 
and NSF DMR 1206728. 


^ E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. 16, 407 
(1961). 

^ M. Oshikawa, Phys. Rev. Lett. 84, 1535 (2000). 

® M. B. Hastings, Phys. Rev. B 69, 104431 (2004). 

^ L. Balents, Nature 464, 199 (2010). 

® M. Rigol and R. R. P. Singh, Phys. Rev. Lett. 98, 207204 
(2007). 

® A. Zorko, S. Nellutla, J. van Tol, L. C. Brunei, F. Bert, 
F. Due, J.-C. Trombe, M. A. de Vries, A. Harrison, and 
P. Mendels, Phys. Rev. Lett. 101, 026405 (2008). 


X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B 83, 
035107 (2011). 

® L. Michel and J. Zak, Phys. Rev. B 59, 5998 (1999). 

® L. Michel and J. Zak, Euro Phys. Lett. 50, 519 (2000). 

L. Michel and J. Zak, Phys. Rep. 341, 377 (2001). 

S. A. Parameswaran, A. M. Turner, D. P. Arovas, and 
A. Vishwanath, Nat. Phys. 9, 299 (2013). 

R. Roy, arXiv: 1212.2944 (2012). 

M. B. Hastings, J. Stat. Mech. 2007, P08024 (2007). 

A. Konig and N. D. Mermin, Proc. Natl. Acad. Sci. 96, 



9 


3502 (1999). 

J. H. Conway and J. P. Rossetti, arXiv:math/0311476 
(2003). 

Our results in this and the next sections are also applicable 
to 2D systems characterized by layer groups, since any 
layer group has a natural correspondence to a space group. 
T. Hahn, ed., International Tables for Crystallography, 5th 
ed., Vol. A: Space-group symmetry (Springer, 2006). 

H. Po, H. Watanabe, M. Zaletel, and A. Vishwanath, 
arXiv: 1506.03816 (2015). 
in preparation. 

R. Lutowski and B. Putrycz, arXiv:1411.7799 (2014). 

J. S. Helton, K. Matan, M. P. Shores, E. A. Nytko, B. M. 
Bartlett, Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, J.-H. 
Chung, D. G. Nocera, and Y. S. Lee, Phys. Rev. Lett. 98, 
107204 (2007). 

Y. Okamoto, M. Nohara, H. Aruga-Katori, and H. Takagi, 
Phys. Rev. Lett. 99, 137207 (2007). 

S. M. Young, S. Zaheer, J. C. Y. Teo, C. L. Kane, E. J. 
Mele, and A. M. Rappe, Phys. Rev. Lett. 108, 140405 
( 2012 ). 

S. M. Young and C. L. Kane, arXiv: 1504.07977 (2015). 

L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

M. P. Zaletel and A. Vishwanath, Phys. Rev. Lett. 114, 
077201 (2015). 

D. Perez-Garci'a, M. M. Wolf, M. Sanz, F. Verstraete, and 
J. 1. Cirac, Phys. Rev. Lett. 100, 167202 (2008). 


Appendix A: Proof of minimal spin-orbit conpled 
filling in the matrix-product-state formalism. 

In the MPS formalism, a short-range correlated trans¬ 
lation invariant wave-function 41 is written as 

X 

E (Ai) 

{a„}=l n 

where the local Hilbert spaces are indexed by a„ and the 
an are “auxiliary” indices to be summed over. For our 
purposes, y can actually be infinite: what’s important 
is that the Schmidt decomposition has a well behaved 
and discrete spectrum, which follows from the area-law 
obeyed by gapped ground states^^. An important ques¬ 
tion is how the global symmetries of 41 manifest them¬ 
selves in the local tensors B. In our case the global U(I) 
and time reversal symmetries factorize into their onsite 
representations. 


n n 

(A2) 

GrK^l[g0^K, 

(A3) 


n 


where n label sites and the g are onsite unitary operators. 
Assuming translation invariance, we will drop the super¬ 
script n. Standard results^^ show that in the absence of 


long-range order, 41 is symmetric if and only if 

Y^geabB^^p = E uL^,B^^,p,Ue p'p, (A4) 

b ot'70' 

YgrabKp = E (A5) 

b 0(',0' 

for some unitary % x % matrices U, and zz S Z is the 
integer filling of the system. 

The transformation laws of Eq. (A5) define the U only 
up to a phase. While the onsite symmetries satisfy 

577 = gror^ (A6) 

since = (—l)^i, the U may realize the symmetry 

group only projectively, 

UrU^ = ojiT, T)U^, to e U(l). (A7) 

In our case this particular relation is not interesting, 
as we can choose the phase of t/,r to set w(T, T) = 1. 
However, when we compare the equalities contained in 
Eqs. (A5) and (A 6 ), we obtain 

g^B = e^^'^UlBU^ = grg*rB = , (A 8 ) 

which implies 

C'^^B = B. (A9) 

Unless v S 2Z, we arrive at a contradiction. 

Appendix B: ‘Modding-out’ and the emergence of 
Spin and Pin^ from a lattice model 

1. Motivation and setup 

In this appendix, we clarify some subtleties in putting 
a system on a compact flat manifold K^/F. To illustrate 
the problem, let us first explore a more explicit form of 
the Hamiltonian on K-^/F obtained by the modding-out 
prescription explained in the main text. 

To that end, let a subspace M C be a fundamental 
domain that contains exactly one of the equivalent points 
of any x G Namely, for any y G there exists a 
unique x G M and 5 G F such that y = g{x). Note that 
g is also unique since F is fixed-point free. Although it 
is not necessary, for simplicity we choose M to be simply 
connected. 

We assume the original Hamiltonian in is a sum of 
local terms H = where Ha^ has a finite support 

Ag around x. Then consider 

Hm=Y 

xGM 

When X in this sum is near the boundary of M, Ha^ 
might contain operators outside of M; in that case they 
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should be understood as the identified operators inside 
M, given by the equivalence relation 

(B2) 

Then the natural identification M ~ Mr = K^/T allows 
us to interpret Hm as the Hamiltonian on Mp. Note that 
the infinitely many difi'erent choices of M gives the same 
Mr as long as one remembers each site in Mr can be 
labeled by the labels of any of the equivalent sites in 

While the choice of g above, as a space group element, 
is unique, for SOC fermions the symmetry group has 
to be extended by the Z 2 fermion parity and g is not 
uniquely defined. As such, a more careful discussion is 
warranted and we tackle it in the last subsection of this 
appendix. 

Putting such complications aside, one still sees there 
are infinitely many valid choices of M. The choice of 
M, however, is merely a gauge choice. More precisely, 
in this appendix we show the following: given any two 
valid choices of fundamental domains M and M', the 
Hamiltonians Hm and Hm' are related by a local unitary 
transformation. 

As an illustrative example, let us consider the following 
Hamiltonian in 2D: 

H= (B3) 

/,mGZ 


where X and Y are spin operators. H is invariant under 
a space group G generated by glide 

: {l,m) !-)■ (l + l, -m), 

GxX(i^ra)G^^ = (B4) 

GxY(^l,m)G^^ = —Yi^l + i-m)' 

and translation Ty, 


Ty : {I, m) !->■ (/, m + 1), 

TyX(l^rn)Ty ^ 

TvY(l^m)Ty ^ = Y^l,m+l)- 


(B5) 


Let r be the subgroup of the full symmetry group gen¬ 
erated by G^^ and Ty'^ with Nx, Ny being odd integers. 
As described, we introduce the equivalence relation : 

~r X(^i m+Ny) ~r (B6) 

Y(l^m) Y(^i^jri+Ny) ~r —Y(^l+N:^-Tn)- (B7) 

For example, M can be chosen as {{l,m) | 1 < Z < 
Nx,l < rn < Ny}. Under identification of the edges 
of M by r, one sees that Mp is actually the Klein bottle. 

As described, the modded-out Hamiltonian is then 
given by 


N,-l Ny 

Nm— ^ ^ ^ ^ ( 1 ) 

1=1 m=l 

Ny-1 

+ X! + ^(iV,,iV„)h"(l,JV„)(B8) 

m—l 


where we have used Eq. (B7) to bring Y(^N^+i,m) into the 
region M and (—= 1 as Nx is odd. 

Staring at Eq. (B8), one sees that: 

1. Around the ‘seam’ 1 = 1, Hm contains couplings 
(boundary conditions) that do not seem to coincide 
with the original Hamiltonian defined on the plane. 

2. The symmetry of Hm is unclear — does Hm still 
have some reminiscent of the G invariance oi HI 

For the first point, we argue that the ‘seam’ of Mp is 
an artifact of choosing a specific fundamental domain. 
We can see that the Hamiltonian (deep) inside M [the 
first sum in Eq. (B8)] is identical to that of the original 
Hamiltonian in By redefining M, any point in can 
be taken far away from the boundary of M. For example, 
for (IojUIq) {IqjItio G If), one may choose 

M = {{l,m) I \l-lo\ < ^,|m-mo| < ^}. (B9) 

Therefore, one can argue that any expectation values of 
local operators or correlation functions around any given 
points in Mp are unchanged from those in the bulk (for 
sufficiently large Nx, Ny). This point can be made clearer 
using the local unitary equivalence we discuss in the next 
subsection. 

The answer for the second point depends on the prop¬ 
erty of F. To understand this, one should ask: does 
[x] = [ 2 ^ always imply [G'(x)] = [G(27)]? It is easy to 
see that this is the case only when F is a normal sub¬ 
group of G. When F is normal, the G invariance of H 
implies the G invariance of Hm (up to local unitaries). 
Even when F is not normal, however, Hm is still invari¬ 
ant under the normalizer of F in the full symmetry group 
(e.g. in the above example, F is generated by Gx only). 
For our purpose, the G invariance of Hm is not necessary 
- the properties that each point of Mp is locally identi¬ 
cal to the bulk and that Hm remains TR symmetric are 
sufficient. 


2. Local unitary equivalence 

Here we explicitly show the asserted local unitary 
equivalence among {Hm} defined with different funda¬ 
mental domains, and discuss its physical consequence. 

To this end, we first define a slightly more formal, but 
nonetheless equivalent, procedure to define the modded- 
out Hamiltonian in Eq. (Bl). This merely serves to intro¬ 
duce notations that make the local unitary equivalence 
manifest. By definition, for any site y in there exists 
X G M and g\j G F, both unique, such that x = g\j(jj). 
In the action of g\j on any operator Oy is repre¬ 
sented by first sending y ^ x and then conjugating Og 
by a local unitary operator s\j. Note that s\j acts non- 
trivially only on site x = g\i(y). After modding out, 
however, the relabeling y x becomes trivial since x, y 
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are merely different labels for the same site. The identi¬ 
fication of operators, as discussed in Eq. (B2), can then 
be defined via the map hm ■ Og i—>■ s^Ogs^. 

can be promoted to a linear map acting on any 
local operators with support in region A C pro¬ 
vided any two distinct sites in A belong to different 
equivalence classes. More concretely, we simply have 

with JJ^j = 

'®M being a local unitary operator. We can now 
rewrite the modded-out Hamiltonian in Eq. (Bl) in a 
more precise manner: 

Hm=J 2 (BIO) 

xGM 

In particular, because the elements of E commute with 
time-reversal T, Hm is time-reversal symmetric. 

Taken at face value, the two Hamiltonians Hm and 
Hm' obtained by modding-out with fundamental do¬ 
mains M ^ M' are not identical. However, the choice 
of the fundamental domain is only a guage choice and 
one can always construct a local unitary transform U 
such that UHmH = Hm'- In fact, one simply sees that 
Hm' = f^M'{HM) and therefore U = Um'- (Note again 
one can freely switch between the equivalent labels for 
the same site.) 

Such local unitary equivalence implies the ground state 
expectation value of any local operators with support 
H.r C Mr is identical to that of the bulk with exponen¬ 
tial accuracy. If Ar corresponds to a region A deep inside 
M, this is evident since the ground state^^ is gapped and 
the restriction of Mm Ur is identical to MU- H cor¬ 
responds to a region A! that crosses the border of M (or 
equivalently the corresponding region is disconnected in 
M), one can nonetheless choose a new fundamental do¬ 
main M' to center A'. In essence, given a ‘measurement 
region’ Ag, one can move any apparent ‘seams’ arbitrar¬ 
ily far away from it using lA (bounded by the size of Mp), 
and hence deduce the measurement outcome is identical 
to the bulk up to exponential corrections. Important to 
the discussion is that U leaves the region Ag unchanged. 
This should be contrasted with a ‘visible’ defect, say an 
impurity, which can be moved only by invoking transla¬ 
tion operator and will therefore shift also the measure¬ 
ment region. 

3. Complications with spinful fermions: Spin and 

Pin 

As explained in the main text, the ‘spin’ degrees of 
freedom complicates the situation since the internal ro¬ 
tations or reflections introduce sign factors that depends 
on the fermion number N. Here we clarify the mathe¬ 
matical structure of Spin or Pin structure on a compact 
flat manifold Mp = M^/T in more detail. 

Before discussing (s)pin, let us first examine the 
Bieberbach group T itself. The first Bieberbach theorem 


states that each 5 S F has a unique decomposition into 
a product of a (possibly fractional) translation T and an 
orthogonal transformation B G 0{D), g = TB. For our 
purposes only a discrete subset G C 0[D) will appear, 
such as reflections and 27: jm rotations. To be fancy, we 
have the short exact sequence 

1-->1, (Bll) 

where 1 / denotes the lattice translations in F and r 
projects onto the point group part of g, i.e., r(TB) = B. 
The map r is the ‘holonomy representation’ of Mp. For 
the Bieberbach manifold Mp, the fundamental group is 
7 r^(Mp) = F; intuitively this is because a straight line 
connecting two points on related by g S F becomes 
a closed curve in Mp. The map r encodes the fact that 
when traversing a path g G 7 r^(Mp), we return to the 
point having been ‘screwed’ by the orthogonal transfor¬ 
mation r(g) G G. 

In the presence of (s)pin degrees of freedom, G is real¬ 
ized as a Z 2 extension encoding the sign structure, which 
we denote by G^. For instance, for spin-1/2, the 27r rota¬ 
tion is G 2 ir = (~l)^j where N denotes the fermion num¬ 
ber. The screw and glide transformations then inherit 
these signs, and a 2 i-screw will square to a translation 
and the fermion parity, Formally, 

1->Z 2 - >G^^^G -^1, (B12) 

where A[(—1)-^] = 1. The most familiar example of 
this sequence is probably the case of G = SO(3) and 
G^ = Spin(3) ~ SU(2). The choice of G^ depends on 
the physical model, but for SOC spin-1/2 fermions G^ 
is a discrete subgroup of the ‘Pin“(M)’ group, the unori¬ 
ented counterpart of Spin(M). The ‘—’ superscript indi¬ 
cates that mirror reflections satisfy = (— 1 )^, since re¬ 
flection is implemented on spins as Rx = icr^ ■ Generally, 
we can consider both Pin'''(M) and Pin~(M) depending 
on whether reflections square to R? = (±1)^. 

To consistently introduce the equivalence relation on 
spinful fermions required to mod-out by F, we need a 
group homomorphism e : F —>■ G^ which encodes how we 
choose the sign that can accompany each reflection and 
rotation. Namely, e has to satisfy 

e(ffi)e(52) = £( 5192 ). (B13) 

such that the following diagram commutes: 

G^ (B14) 

e / 

/ A 

F^^^^G 

This is a ‘lift’ of the holonomy representation to G^. De¬ 
pending on F, it may or may not be possible to find such 
an e. In fact. Ref. 20 shows that Eq. (B14) is precisely 
the condition for the existence of a Spin / Pin structure 
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on Mr- All orientable 3-manifolds admit spin structure 
and all nonorientable 3-manifolds admit Pin” structure. 
Hence, all 10 Bieberbach 3-manifolds (6 orientable and 4 
nonorientable) admit Spin or Pin” structure, and such e 
always exists. 

As an illuminating example, let us discuss the space 
group No. 29. This group is generated by three symmetry 
elements: 


Gx ■■ (x, y,z) ^ (x + ^, -y, z), 

(B15) 

Ty : {x,y,z) 1 -^ {x,y + l,z), 

(B16) 

Sz ■■ {x, y, z) ^ {-X, -y, z + ^) 

(B17) 

with relations 


^2 _ rp q2 _ rp Q _ 'T' Q 

- -^215 '^z - “^2’ X>^Z^X ^ 

(B18) 

where Tx^y^z is translation by one unit. Since 
IG 2 X {I is the inversion) and r{Sz) = G 22 , we 
to signs ^x,y,z = ±1: 

r(Gx) = 
have e up 

e(Gx) = 

(B19) 


(B20) 

e(Sz) = ^ze^^^ =Czaz. 

(B21) 


For e to be homomorphic, e has to keep relations in 
Eq. (B18). It is easy to see that once we fix e{Tx) = 
e{Tz) = —1, all relations are satisfied without specifying 
^x,y,z = ±1- We thus have 2^ = 8 Pin” structures. 

Curiously, this group, No. 29, is the single exception 
that does not admit Pin^ structure among all Bieberbach 
groups. For Pin'*" we have to use e(Ga;) = so that 
e(G^) = -1-1, which also fixes e(Tx) = -1-1. Then, the third 
relation in Eq. (B18) is violated. Therefore, we cannot 
put particles obeying Pin'*' on the manifold ‘1st Amphidi- 
cosm’ obtained by modding-out by F = No. 29. This, 
however, does not concern us here since we are interested 
in spin-1/2 electrons. 

Note that because we consider models with U(l) sym¬ 
metry, it seems that we could redefine the reflections by 
an arbitrary U(l) phase arbitrary phase 6, —)■ Rx- 


Choosing then 9 = tt/2 erases the distinction between 
Pin'*’ and Pin”. Formally this amounts to considering 
U(l), rather than Z 2 , extensions of the space group, lead¬ 
ing to Spin^ and Pine structures. However, time-reversal 
narrows down this U(l) ambiguity since the modding-out 
process must preserve TR symmetry. The space group 
elements g € T and T commute; the lift e must as well, 
^{ 9 )'^ = '^'^( 9 )^ otherwise the ‘seams’ of the region will 
break TR. For S' = 1/2 models, this amounts to requiring 

<^y<^i9y<^v = 49)- (B22) 

So TR invariance reduces the U(l) ambiguity back down 
to Z 2 , fixing either Pin^. For example, e{Gx) = ^xi<^x for 
Pin” satisfies this, while e{Gx) = ^x^x for Pin"'' does not. 
In general, the relation R^{TRT~^R~^) = (±1)^ is in¬ 
variant under U(l) redefinitions of R (since T conjugates 
phases) and determines which of Pin^ to use. 

Appendix C: Group-subgroup relations and 
elementary NS space groups 


As explained in the main text, all 157 NS space groups 
contains at least one of the elementary NS space groups 
as t-subgroup. Figure 4 below illustrates which space 
group contains which elementary space groups, gener¬ 
ated based on Ref. 17. Elementary NS space groups are 
grouped into 10 Bieberbach manifolds^®. 

Using this group-subgroup relation, one can deduce the 
constraint on filling v for all non-elementary space groups 
G based on the result for NS elementary space groups es¬ 
tablished in the main text. It can in general be expressed 
as v = mgn, where 


mg = max 


\G--T\ 

fo<6 IS : Fol 


(Cl) 


with Fq being a Bieberbach (fixed-point free) subgroup 
of G and T being the lattice translation subgroup of G- 
(Any Bieberbach group is ‘elementary’ in our notation, 
but some elementary space groups are not fixed-point 
free.) We list the results in Table H. 
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TABLE II. The filling that may realize sym-SRE phases for nonsymmorphic space groups. Those space groups not listed here 
are symmorphic and hence v = 2n; n in this table is a positive integer 1, 2, 3, • • • . 


No. 

No. 

No. 

No. 

No. 

No. 

No. 

4 

4n 

39 

4n 

66 

4n 

100 

4n 

129 

4n 

165 

4n 

201 

4n 

7 

4n 

40 

4n 

67 

4n 

101 

4n 

130 

8n 

167 

4n 

203 

4n 

9 

4n 

41 

4n 

68 

4n 

102 

4n 

131 

4n 

169 

12n 

205 

8n 

11 

4n 

43 

4n 

70 

4n 

103 

4n 

132 

4n 

170 

12n 

206 

4n^ 

13 

4n 

45 

4n 

72 

4n 

104 

4n 

133 

4n^ 

171 

6n 

208 

4n 

U 

4n 

46 

4n 

73 

4nt" 

105 

4n 

134 

4n 

172 

6n 

210 

4n 

15 

4n 

48 

4n 

74 

4n 

106 

4n^ 

135 

4n^ 

173 

4n 

212 

8n 

17 

4n 

49 

4n 

76 

8n 

108 

4n 

136 

4n 

176 

4n 

213 

8n 

18 

4n 

50 

4n 

77 

4n 

109 

4n 

137 

4n 

178 

12n 

214 

4n 

19 

8n 

51 

4n 

78 

8n 

110 

4n^ 

138 

8n 

179 

12n 

218 

4n 

20 

4n 

52 

8n 

80 

4n 

112 

4n 

140 

4n 

180 

6n 

219 

4n 

24 

4n 

53 

4n 

84 

4n 

113 

4n 

141 

4n 

181 

6n 

220 

4n'^ 

26 

4n 

54 

8n 

85 

4n 

114 

4n 

142 

4n^ 

182 

4n 

222 

4n 

27 

4n 

55 

4n 

86 

4n 

116 

4n 

144 

6n 

184 

4n 

223 

4n 

28 

4n 

56 

8n 

88 

4n 

117 

4n 

145 

6n 

185 

4n 

224 

4n 

29 

8n 

57 

8n 

90 

4n 

118 

4n 

151 

6n 

186 

4n 

226 

4n 

30 

4n 

58 

4n 

91 

8n 

120 

4n 

152 

6n 

188 

4n 

227 

4n 

31 

4n 

59 

4n 

92 

8n 

122 

4n 

153 

6n 

190 

4n 

228 

4n^ 

32 

4n 

60 

8n 

93 

4n 

124 

4n 

154 

6n 

192 

4n 

230 

4n^ 

33 

8n 

61 

8n 

94 

4n 

125 

4n 

158 

4n 

193 

4n 



34 

4n 

62 

8n 

95 

8n 

126 

4n 

159 

4n 

194 

4n 



36 

4n 

63 

4n 

96 

8n 

127 

4n 

161 

4n 

198 

8 n 



37 

4n 

64 

4n 

98 

4n 

128 

4n 

163 

4n 

199 

4n 




^ t- = 8n — 4 is prohibited for the noninteracting case^®. There is no known interacting model of sym-SRE at these filling either, 
y = 4 is prohibited for the noninteracting case^®. There is no known interacting model of sym-SRE at this filling either. 
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Classification of 
all 230 space groups 

Symmorphic 

r 


Name of flat manifold 


ITC No. (simbol) of elementary group F / condition on v. 

List of space groups whose maximal elementary f-subgroup is L. 

The bold italic indicates that the most symmetric Wyckoff position of 
the space group corresponds to a lattice with more than v^J2 sites in each 
unit cell, implying that cannot be realized by any atomic insulator. 


Toms 


No. 1 (PI) G 2Z 


1, 2, 3, 5, 6 , 8 , 10, 12, 16, 21, 22, 23, 25, 35, 38, 42, 44, 47, 65, 69, 71, 75, 79, 81, 82, 83, 87, 89, 97, 99, 107, 
111, 115, 119, 121, 123, 139, 143, 146, 147, 148, 149, 150, 155, 156, 157, 160, 162, 164, 166, 168, 174, 175, 
177, 183, 187, 189, 191, 195, 196, 197, 200, 202, 204, 207, 209, 211, 215, 216, 217, 221, 225, 229. 


Non- symmorphic 

^-Dicosm-^ 

No. 4 (P2j) IV e 4Z 

4, 11, 14, 17, 18,20, 26,31, 
36,51,53,55,58, 59, 63,64, 

90, 94, 113, 114, 127, 128, 

129, 135, 136, 137, 173, 176, 
182, 185, 186, 193, 194. 

V___ J 

. -Didicosm-. 

pNo. 19 (P2,2j2j) iv G 8 Z ^ 

19, 61, 62, 92, 96, 198, 
205,212,213. 


p No. 24 (/2j2,2,) tv G 4Z ^ 

24, 75, 74, 98, 122, 141, 
142,199, 206, 210, 214, 
220, 227, 228, 230. 

V 


y 


-1st Amphicosm- 

No. 7 (Pc) IV G 4Z 

7, 13, 14, 26, 27, 28,30,31,32, 
34,39,41,48,49, 50,51,53, 
55, 58, 59, 64, 67, 68 , 85, 86 , 
100, 101, 102, 103, 104, 106, 
116, 117, 118, 124, 125, 126, 
127, 128, 129, 132, 133, 134, 
135, 136, 137, 201,222, 224. 


.v- -Tetracosm- s 

r No. 76, 78 (P4,, P 43 ) IV G 8 Z n 

76, 78, 91, 92, 95, 
96,212,213. 


-No. 77 (P 43 ) IV G 4Z ^ 

77, 84, 86 , 93, 94, 101, 102, 
105,706, 131, 132, 133, 134, 
135, 136, 137, 208, 223, 224. 


-No. 80 (/4j) tv G 4Z- 

80, 88,98, 109, 110, 141, 

142, 210, 214, 221, 228, 230. 

V y 


X-2nd Amphicosm- 

No. 9(Cc) tvG4Z 

9, 15,36,37, 40,41,43,45, 

46, 63, 64, 66 , 68 , 70, 72, 75, 

74, 88 , 103, 104, 105, 106, 

108, 109, no, 112 , 114, 120 , 
122, 124, 126, 128, 131,755, 
755, 137, 140, 141,742, 158, 
159, 161, 163, 165, 167, 184, 
185, 186, 188, 190, 192, 193, 
194, 203,206,218,219,220, 
222, 223, 226, 227, 228, 230. 

\ ___ J 


Tricosm 


No. 144, 145 (P 3 ,,P 33 ) iv G 6 Z 


144, 145, 151, 152, 153, 
154, 171, 172, 180, 181. 


y 


-Hexacosm-^ 

No.169,170 (P6j,P63) v G 12Zn 
169, 170, 178, 179. 


r No. 171,172 {P6^,P6^ IV G 6 Z^ 
171, 172, 180, 181. 

-No. 173 (P 63 ) IV G 4Z- 

173, 176, 182, 185, 186, 

193, 194. 

V_^_ J 


— 1 st Amphidicosm — 

No. 29 (Pcfl2j) tv G 8 Z 
^ 29, 54, 57, 60, 61, 205. ^ 


- 2nd Amphidicosm - 

No. 33 (Pfl2j) IV G 8 Z 
33, 52, 56, 60, 62, 130, 138. 


FIG. 4. Classification of all 230 space groups based on its associated Bieberbach manifolds. 


















































































